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The corrections of order (Za) 6 mi/ni 2 and (Za) 7 from one-loop two-photon exchange diagrams to 
the energy spectra of the hydrogenic atoms are calculated with the help of the Taylor expansion of 
corresponding integrands. The method of averaging the quasipotential over the wave functions in the 
d-dimensional coordinate space is formulated. The numerical values of the obtained contributions 
to the fine structure of muonium, hydrogen and positronium are presented. 


I. INTRODUCTION 


The study of perturbative effects of higher order in a ( a is the fine structure constant) in the energy spectra 
of the positronium and muonium is of great practical importance for testing the bound state theory in quantum 
electrodynamics. In the past few years considerable progress has been made towards that goal [jl], First of all it is 
determined by an increase of the accuracy in the measurement of the fine and hyperhne structure in these leptonic 
systems. Thus, the experimental error for the muonium ground state hyperfine splitting was reduced by a factor of 
three 101: 


A^hfs(Mu) = 4 463 302 765(53) Hz 


( 1 ) 


Even larger increase of the precision was reached in the measurement of the interval IS 4- 2S in the muonium on the 
basis of the Doppler-free two-photon spectroscopy ||: 

(2 3 Si 4 1 3 Si) = 2 455 528 941.0(9.8) MHz. (2) 

The value of the same interval in positronium obtained some years ago is the following [Q : 


Ai/pg P (2 3 Si 4- l 3 Si) = 1 233 607 216.4(3.2) MHz. 


( 3 ) 


One of the most precise experiment in the spectroscopy of simple atomic systems is the measurement of the hydrogen 
gross structure IS 4- 2S ||: 


Ai/h( 2S - IS) = 2 466 061 413 187.34(84) kHz 


( 4 ) 


On the other hand this progress results from the development of the computer methods for the calculation of Feyn¬ 
man amplitudes and the appearance of nonrelativistic quantum electrodynamics (NRQED) [^,^] for the calculation 
of the bound state energy spectra at the nonrelativistic scale. Due to performed calculations there were obtained some 
new contributions to the energy levels of hydrogen-like systems: corrections of order (Zcc) 6 mi/m 2 to the muonium fine 
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structure |pl. [TTl| . |TT| . |Tl^ . |TTs|JT~i| l , corrections of order ma 6 to the positronium fine and hyperfine structure [ 0 , |l6l0 . 0 , |l9| , 
logarithmic contributions 0(ma 7 ln 2 a) to the positronium spectrum |2(j| . ^T|] and also the logarithmic corrections of 
order 0(a' In a) to the muonium and the positronium hyperfine splittings ]22|j2^| . Calculation of loop integrals which 
represent the corrections to the Coulomb potential by the perturbation theory is complicated problem because of 
the presence of the essentially different energy scales determined the behaviour of the integrands. The following four 
scales can be introduced in terms of the loop energy p° and the loop three momentum p ^4|,^5|] : 


hard momentum region : 

: P° ' 

- /.(, p ~ p, 

(5) 

potential region : p° ~ 

pa 2 

, p ~ pa, 

(6) 

soft momentum region : p 


pa, |p| ~ pa, 

(7) 

ultrasoft momentum region : 

P°^ 

- pa 2 , p ~ pa 2 , 

(8) 


where the mass parameter p is determined by the masses of the particles mi, m 2 . There are two approaches to the 
calculation of the definite order contribution on a. the first one is related with the explicit extraction of the small 
parameter in the correspondi ng QE D Feynman amplitudes which than can be expanded in Taylor series in one of 
the energy regions (5)-(8) 0 . The second approach realized in NRQED consists in the formulation of such 

procedure already at the Lagrangian level . In present work which is the sequel of [|l4| we study the contributions 

of order (Za) 6 and (Za) 7 to the spectra of muonium and positronium from one-loop Feynman amplitudes using the first 
approach. The method of dimensional regularization and Taylor series allows to construct the interaction operator 
of the particles in the momentum and the coordinate representations. Numerically the order factors are equal to 
(Za) 6 m 2 /m /J =0.0902 MHz and m e a'=0.136 MHz so the calculation of the mentioned contributions is very actual 
especially for their confrontation with experimental data (1) and (4). 


II. CONTRIBUTIONS OF ORDER (Za) 6 Mi/M 2 TO THE MUONIUM FINE STRUCTURE. 


Simple counting of the Za powers in the two-photon exchange diagrams in Fig.l shows, that they can contribute 
in the order (Za) 2 • (Za) 3 = (Za) 5 . Here the first factor is connected with exchanged photons and the second one 
is determined by the wave functions of the bound state. At the same time the availability of the relative motion 
momentum of the particles in the initial state p (|p| ~ pa) and in the final state q (|q| ~ pa) in these 2y amplitudes 
lead to the appearance of the higher order contributions on a. First of all we formulate here the approach to the 
calculation of definite order corrections and then calculate the contributions containing additional degree Za and 
the ratio of the electron to muon masses m 1 /m 2 = 0.004836. This approach can be used also for the positronium 
spectrum if we don’t employ some simplifications in the 2y Feynman amplitudes connected with small value nij /m 2 . 
The general structure of the integrals describing the quasipotential remains unchanged in this case. 

Our calculations are based on a local quasipotential equation of the Schroedinger type 0]: 

(9) 


where b 2 = E 2 — mf=E| — nij, pr = E]E 2 /M is the relativistic reduced mass, M = Ej + E 2 is the mass of the bound 
state, iPm(p) is the quasipotential wave function. For the initial approximation of the quasipotential V(p, q, M) for 
the bound system we take the ordinary Coulomb potential. 

To construct the interaction operator of the particles corresponding to 2y amplitudes we used also the projection 
operators 7rs=o and 7rs=i on the states with total spin S=0 and 1 in the system (e _ p + ): 


P + M „ P + M, 
~ 2 V / 2M 75 ’ 7FS=1 ~ 2\/2M 6 


where is the polarization vector for the state 3 Si, e = 


( 10 ) 


2 





Pi = ( E 1 ,p ) 


(Ei + fc 0 , k) 


qi = (E 1 ,q) 




P 2 = (E 2 , -p) (E 2 - ko, -k) q 2 = (E 2 ,~ q) 

FIG. 1. Two-photon exchange Feynman diagrams. P = pi + P 2 is the total momentum of the two-particle bound state, p,q 
are the relative four-momenta of particles in the initial and final states. 


The contribution of diagram (a) Fig.l takes the form (we consider further only triplet states 3 Si): 

( a) i(Za) 2 r f 1 (k,m 1 ,m 2 )d 4 k 

v 27 lP,qj ^ 2 J (k _ p) 2 (k _ q) 2 Dl(k)D2( _ k) > 


( 11 ) 


fi(k, m 1; m 2 ) = m 2 (4mi + 2k°) + 2mik° - 2k° 2 + ^k 2 . 

o 

The denominator of the electron propagator D;i(k)= k 2 + 2m 1 k° + b 2 . But corresponding muon denominator can be 
simplified taking into account the small parameter mi/m 2 and neglecting by the muon kinetic energy in the interme¬ 
diate state: D 2 (—k) ss 2m 2 (—k° + ie). The crossed two-photon diagram (b) Fig.l gives the similar contribution: 


( b ), , _ i(Za) 2 I" _ f 2 (k, mi, m 2 )d 4 k _ 

2t p,q 7T 2 J (k-p) 2 (k-q) 2 Di(k)D 2 (-k + p + q)’ 


( 12 ) 


9 10 

f 2 (k,mi,m 2 ) = m 2 (4m 1 + 2k°) + 2mik° + 6k°" + —-(pk + qk - k 2 ). 

o 

Direct two-photon amplitude (a) Fig.l contains the contribution of lower order on a (~ a 3 ) which is cancelled by 
the iteration part of the quasipotential Vi 7 x G f x Vi 7 |l4|]. To extract such term we transform the product of the 
electron and muon denominators as follows ]25[ |: 

1 —27ri<5(k°) 1 

Di(k)D 2 (-k) “ —2E(k 2 — b 2 ) “ 2E 


(k° + iO)Di (k) (-k° + iO)D 2 (-k) 


(13) 


Omitting the first term in the right part of (12) corresponding to free two-particle propagator we can present the sum 
of potentials V 2 ^ +V 2 *^ in the form: 


V 27 (p,q) = (p, q) + (p, q) 


i(Za) 2 

6m 2 7r 2 


d 4 k 

(k — p) 2 (k — q) 2 X 


(14) 


12m 2 + 18mik° + 24k° 2 - 12k 2 + 10k(p + q) 6mi + 3k° 
(k° + iO)Di(k) (-k° + iO) 2 


Consider first of all the necessary order contribution of (14) in the potential region (6). In this case the integration over 
k° is determined by the residue in the pole of the electron propagator, which can be presented as Taylor expansion: 


1 

k° 2 — k 2 + 2mik° + b 2 


E 


n—0 


(-k° 2 ) n 

(—k 2 + 2mik° + b 2 ) n+1 


(15) 
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The expansion of photon propagators in the potential region looks as follows: 


1 


= £■ 


(-k° ) n 


k°"-(k-p) 2 ^ Q [-(k-p) 2 


)2]n+l ‘ 

Multiplying the expansions (15) and (16) we observe the next structure of (14): 

1 k ° 2 k ° 4 


V£°‘(p,q) ~ J d 4 k 


-(k-p) 2 (k-p) 4 (k-p) 6 


1 


k° 


k ° 4 


— (k — q) 2 (k-q) 4 (k - q) 6 


1 


(k° + iO) 


( 16 ) 


(17) 


1 _k^_k°^_ 

—k 2 + 2 mik° + b 2 (—k 2 + 2 mik° + b 2 ) 2 (—k 2 + 2 mik° + b 2 ) 3 

The integral (17) contains corrections of different order on a which is determined by the corresponding factors after 
integration over k° in the pole of electron propagator. Existing contribution of order (Z a) 6 in Eq. (17) is determined 
by the following typical integral: 


J(p,q) 


dk _ (k 2 ) a _ 

(27r) 3 (k 2 — 2kp + p 2 )^(k 2 — 2kq + q 2 ) 7 ’ 


(18) 


where the degrees a, /?, 7 can be different, but the sum 2a + 3 — 2/3 — 2y = 1. The contribution to the energy spectrum 
of this system can be obtained after averaging J(p, q) over Coulomb wave functions. The characteristic integrals (18) 
are divergent just as in ultraviolet region so also in infrared one. The reason of these divergences lies in the used 
Taylor expansions of the integral function in Eq. (14) which we made to extract the necessary order contribution 
on a. The initial integral (14) is finite so all divergences in the intermediate expressions of the type (18) must be 
cancelled in the sum. There are two ways for the calculation of the average value < J(p, q) >: 

1. The integration (18) over k can be done by means of the dimensional regularization. Extracting the finite part 
of the result we can construct the Fourier transform and averaging it on Coulomb wave functions in the r-space. 

2. Carrying out the Fourier transform of Eq. (18) directly we can then average it in coordinate space. 

We used both approaches for the calculation of the contributions (Za) 6 mi/m 2 in the muonium energy spectrum. The 
transition to the d-dinrensional space allows to regularize the expression of the quasipotential both in the relativistic 
region of the intermediate momenta and in nonrelativistic one. The results of the integrations (18) for different 
functions are presented in the Table I. 

The region of the soft momenta (7) gives the necessary order contributions to the muonium S-states. In this case the 
value of the integral (14) is determined by the residue in the photon poles. The expansion of the electron propagator 
in the region (7) takes the form: 


k ° 2 - k 2 + 2 mi k° + b 2 ^ ( 2 mik °) n + 1 


_ g (k ° 2 - k 2 + b 2 y 


(19) 


Using it we have the following standard integrals in Eq. (14): 


I so ft(p,q) ~ J! 


cl 4 k 


1 (k ° 2 -k 2 +b 2 ) E 


n=0- 


(k — p) 2 (k — q ) 2 (k° + iO) (2mik °) n+1 


( 20 ) 


Since the function in the numerator of Eq. (14) contains the terms of order 1, a, a 2 , the sum (20) gives the 
contributions of order 0(a 6 ) when n=l, 2, 3. The subsequent calculation of such integrals was carried out as in the 
potential region ( 6 ). 
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III. MATRIX ELEMENTS OF THE OPERATOR R~ v IN D-DIMENSIONAL COORDINATE SPACE. 


In many cases the averaging of the two-photon interaction quasipotential on Coulomb wave functions can be 
conveniently done in coordinate representation. Indeed the calculation of the integrals J(p,q) shows that the matrix 
elements < J(p,q) > in three dimensional momentum space are divergent. The reason for these divergences is the 
same: Taylor expansions of the integrands. We can implement again the dimensional regularization as for the 
integration on loop momentum k. But then there is need to know the two-particle wave function for the Coulomb 
potential in d-dimensional momentum space J29|,[30]] . Other approach for the calculation of the average value < V 2 7 > is 
connected with the construction of Fourier transform for the potential (14) in d-dimensional coordinate representation 
and subsequent calculation of the matrix elements on Coulomb wave functions in r-space. The Fourier transform of 
the power potential V(p) = l/p m is determined by the following expression j3ij : 


d-3 


V(r) = 


(2tt) 


d- 


_ r(^±i)r(i/2) 


k=l 


r( 2/ 


sin 


d—2 


9e ipr 


COB ®d 0 p d-1-r 


dp = 


F(d^n) 


'0 JO 


2m 7r d/2 r d-m r (m)' 


( 21 ) 


So in three dimensional coordinate space d=3 we have the degree potentials which lead to the divergent matrix 
elements for the case of S-states. To extract such singular terms of the kind 1/e (d = 3 + 2e) we formulate auxiliary 
Coulomb task in d-dimensional coordinate space. Let consider the ordinary Coulomb potential in d-dimensional space 
01 : 


H =T +Vc « = f + i 


■ 2 , (d — l)(d — 3) 


L + 


-, a = p,Za, r = . 

r \ 


E : 

i=l 


( 22 ) 


where p r is radial momentum operator, fj, is nonrelativistic reduced mass. The proper functions of the angular part 
L 2 of the Laplace operator A in d-dimensional space are the homogeneous harmonic polynomials of the degree 1 on 
sphere S d_1 with proper values A = 1(1 + d — 2) j32|. The exact solution of the d-dimensional Coulomb problem in 
the case of discrete spectrum was obtained in Ref. jdij. In particular it was shown that Coulomb energy levels are 
the following: 


e n 


1 


2 (n + 


^) 2 ’ 


(23) 


where n=l, 2, 3, ... is the principal quantum number (see Appendix A). The Coulomb Green function in d-dimensional 
case is important for the calculations of the corrections in the higher orders of the perturbation theory. The expression 
for it was obtained in Ref. f33| . The reqqurence relation for the matrix elements of the power operators r~ v are very 
important for finding the average values of the quasipotential V 2 7 . To obtain this relation in d-dimensional space 
consider the operator of the radial momentum p r : 


Pr = 


1 f r i 


/ ri \ ri 1 
\7’ Pi } = T Pi ~~ 


2 lr ) r 
It satisfies to the following commutation relations: 

[r, p r ] = i • I, 

From the expressions (24), (25) it appears that: 


( d — 1 ) (d d -1 

2 r V dr 2 r 


^7’Pr 


pZZ+1 


H, 


1 


v- 1 


2 i \ . A 0 a 2 (d—l)(d —3) 

- Pr ), [H,i Pr ] = -^--, A 0 = L 2 + --^- L. 


Another important relation can be obtained from (24)-(26) 




1 


v.U +1 


,H 


2v v + 1 

-H- — 

•piz+3 




(d-l)(d-3) 


(24) 


(25) 


(26) 


(27) 


{2v + l)a v{v + l)(v + 2 ) 1 
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So the reqqurence relation for the matrix elements ~ can be derived after averaging (27) on the wave functions 'tp n \ m 
satisfying the following properties: 


H^nlm — E n — ci € n — 


2 (n - 


d-3\2 ’ 


L 2 Vvim = 1(1 + d - 2)i/j nhu . 


(28) 


The matrix element of the left part (27) on ip n i m is equal to 0. So we obtain the necessary reqqurence relation in the 
form: 


0 = 2uE n < — >-(„ + !) 


1(1 +d-2) + 


(d — l)(d — 3) 


< 




> 


(29) 


-|-a(2i/ -|- 1) <C 




> 


v(v+l)(v + 2) 1 


< 


-v+3 


> 


In the case of S-states 1=0 and the obtained reqqurence relation (29) can be rewritten as follows: 

1 8 E n ^ ^ 1 

< ^7+3 >“ (i/ + l)[+(r/ + 2) — (d — l)(d — 3)] < > - 


(30) 


4a(2u + 1 ) 1 

(v + l)[v(v + 2) — (d — l)(d — 3)] < r "+ 2 > 

This formula allows to express the matrix elements of high negative powers of the radius p which are singular in d=3 
space through the matrix elements of lower powers r extracting explicitly the singular factors - (d = 3 + 2 e). As an 
example for the application of (30) we consider the calculation of the typical integral in (14) using two independent 
approaches: 


i(p,q) 


dk ( 47 ra) 2 (k 2 — b 2 ) 
( 2 ^)d (k-p) 2 (k-q) 2 2 / 


(31) 


where the factors ( 47 ro :) 2 and 1/2 fj, were introduced to simplify the intermediate expressions. Indeed the Coulomb 
wave function satisfies to the Eq. (9) with Coulomb potential [^5|. Using it we can transform (31) to the sum of the 
different matrix elements: 


< i(p,q) >=< 


/k 2 -b 2 

V 2 /r 


>=< (k 6 + 3k 4 W 2 + 3k 2 W 4 + W 6 ) 


8 /x 3 


>, 


W 2 = -b 2 . 


Some of them can be calculated directly in d=3 space: 


< k 2 >=< 2/r(E„ - V c ) >= 

n z 

< k 4 >=< (E n - V c )(E n - V c )4 M 2 >= 


< W 6 


Qq ) 6 

n 6 



But the different situation arises when we average k 6 : 

< >=< (En - Vc)^(E n - Vc) >=< E 2 (E n - Vc) > - 

8 fi* 2 n 


(32) 


(33) 

(34) 


(35) 


W 2 a 2 

-2 < E„(E„ - Vc)V c > - < —\ 2 C > — < V^ > + < > • 

There are two divergent matrix elements in three dimensional space: ~< p > and < p >. Using our reqqurence 
relation (30) we observe that the singular parts of (35) at are cancelled in the limit d —> 3. We have the following 
finite result for it: 


6 






















( 36 ) 


k 6 

Q « 


5 

7 

8 ' 

8n 6 

3n 5 

’ 3^ 


Summing the contributions (33)-(34), (36) in the expression (32) we find the value of the initial integral: 


< i(p,q) >= /t 3 « 6 


3n 5 3n 3 


(37) 


Consider the alternative method for the calculation of (31). We can carry out direct integration in (31) by means 
of dimensional regularization. Then going to the d=3 we obtain the following contribution to the quasipotential in 
momentum space: 


i(p,q) 


(47ra) 2 p 2 — b 2 
2 fi 8t 


t 

16 


,t = (p- q). 


(38) 


The Fourier transform of (38) allows to derive the corresponding potential in the coordinate representation: 


I(r) 


2M [ 


1 1 


(P + 


(39) 


Average values of the terms in this relation can be found using the equation for the wave function and reqqurence 
relation (30): 


< 


(p 2 — b 2 ) a 2 
2/r r 2 


3 „ ,6 


>=< -r- >= /T a 


6 6 


en 0 m* n J 


< 


2fir A 


>= 


/i 3 a 6 


en J 


20 

Tn 3 


12 

— r + 


3n 5 


(40) 

(41) 


where the contributions ~ 1/e in two last divergent matrix elements are distinguished only by the sign. As a result 
the sum of (40) and (41) is finite and coincides with (37). Both methods for the calculation of the integrals similar to 
I(p, q) were applied in this work to obtain the contribution of order (Z«) 6 mi/ni 2 in the muonium spectrum. Total 
sum of all contributions in the potential region (6) and soft region (7) to the energy levels is equal 


< v 2 , >= 

m2 


37 1 7 

9n 3 6n 4 9n 5 


(42) 


TABLE I. Results of the integrations (18), t = (p — q). 


1 a 

0 

l 

2 

II 

C" 

i—l 

II 

1 

p 2 +q 2 -t 2 


8t 

16£ 


/3=1,7 = 2 


p 2 —q 2 +t 2 

p 4 -3(q 2 -t 2 ) 2 +2p 2 (q 2 +t 2 ) 


16t 3 

32t 3 

co 

II 

C- 

1—1 

II 



3p 4 +3(q 2 -t 2 ) 2 +p 2 (-6q 2 +2t 2 ) 



64 1 5 

0 = 2,7 = 2 

— 

— 

t 4 + 2t 2 (p 2 +q 2 ) — 3(p 2 — q 2 ) 2 

32 1 5 
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IV. CORRECTIONS OF ORDER ( Za ) 7 TO S ENERGY LEVELS. 


As was pointed out above the main contribution of two-photon amplitudes to the energy spectrum is of order 0(a 5 ). 
These amplitudes give also the contribution of order (Za) 7 in the hard momentum region. To obtain such corrections 
to the S-levels of the bound state we average 2y amplitudes over the particle spin. As a result the 2y quasipotentials 
(a) and (b) Fig.l can be presented in the form: 


(a). . = i(Za) 2 f gi(k,mi,m 2 )d 4 k 

v 27 lP,qj J (k _ p) 2 (k _ q) 2 Dl(k)D2( _ k) > 


(43) 


gi(k,mi,m 2 ) = 16mim 2 — 8mik° + 8m 2 k° — 4W 2 ( —- H-- ) + 

V m 2 mi ' 


V 27 (p,q) = 


+8k u W 2 ( —-— ) - 16k u " + 8k 


i(Za) 2 


mi m 2 / 

g 2 (k,mi,m 2 )d 4 k 


^ 2 J (k — p) 2 (k — q) 2 Di(k)D 2 (—k + p + q) ’ 

g 2 (k,mi,m 2 ) = 16mim 2 + 8mik° + 8m 2 k° - 4W 2 ( — + — ) - 

\m 2 mi/ 


(44) 


—8k°W 2 f —-— ) + 16k° 2 - 8k 2 + 8kp + 8kq, 

V mi m 2 / 

Di, 2 (k) =k 2 + 2Ei, 2 k°-W 2 . (45) 

Taking into account that k° ~ fi, |k| ~ n the expressions of the electron and photon propagators are the following: 


= E 


(W : 


2 \n 


D i, 2 (k) ^ (k 2 + 2Ei, 2 k°) n +i' 


(46) 


1 1 P 2 - 2pk 4(kp) 2 


(k-p) 2 k 2 (k 2 ) 2 


(k 


213 


(47) 


where we kept only the terms of the necessary order in the last expression. Substituting (45) and (47) to (43), (44) 
we revealed that the hard part of the quasipotential V 27 is determined by the set of the characteristic integrals: 

f d 4 k S(k) 

K(p, q) — / 77777TTTvTnTT - ’ (^®) 


16i7T 2 (k 2 )«(k 2 + 2k°Ei)/ 3 (k 2 ± 2k°E 2 )^ ’ 


TABLE II 


. The results of some basic integrals (48): / ^ (t2)a(t2+2t .^ (t 2 _ 2t » E2) 


S(k) 

1 

k° 

k 2 

k° =* 

II 

c" 

T—i 

II 

n. 

of 

II 

S 

e 1 e 2 -e2-e% 

ft 

1 

ft 

to 

1 

— 

96 E'(E% 

32 E'(E% 

2>“2E\ E 2 

« = 2,/3 = 2 , 7 == 1 



— 

e2-e%-e 1 e 2 



128E'(E%(E 1 +E 2 ) 

a = 3,/?=l , 7 = l 

— 

(El+El){E 1 -E 2 ) 

E^ — E 1 E 2 +E 2 

E^—Ei E 2 +E 2 

192 E*E* 

144 E'(E% 

64 E'(E% 




























where the function S(k) contains the powers of k° or (kp). The dimensional regularization was also used for the 
calculation of the integrals (48). The corresponding results are presented in the Table II where we wrote only the 
finite parts of the obtained expressions. 

Summing all contributions in (43), (44) after integration (48) we find that the spin-independent part of the two- 
photon quasipotential has the following structure: (Za 2 p 2 F(mi, m 2 ). The contribution of this term to the energy 
spectrum can be obtained by the use of the dimensional regularization and the equation (9) : 

< p 2 >= J ^^-p 2 ^(p) = b 2 |V>M(0)| 2 . (49) 

So two-photon diagrams (a) and (b) Fig.l yield the following correction of order (Za) 7 for the positronium spectrum 
in the hard momentum region: 


AB„(Ps) = 


( Zq 0 7 

7 m 5 


mi 


1129 
720 ' 


(50) 


In the case of the muonium we can make additional expansion on the powers mi/m 2 . Preserving only linear terms 
on mi/m 2 we obtain: 


AB„(Mu) = -<^m 1 (| + 20iHi). (51) 

7rn° \ 5 m 2 / 

We took into account here only power corrections (Za) 7 . The logarithmic terms of the kind h^mi^/A) (A is the 
parameter of the dimensional regularization) were omitted systematically while they are present in the intermediate 
expressions (48). 

In the soft momentum region the general structure of the quasipotential which gives the corrections of order (Za) 7 
is similar to (20). In the leading order in mi/ni 2 the muonium quasipotential has the following form: 


V 


soft 

27 


(p,q) 


i(Za) 2 

47r 2 mfm2 


d 4 k [ (k 2 - W 2 ) 4 

(k-p) 2 (k-q) 2 k 2 \ k 4 + 


(52) 


k 2 - W 2 ) 3 (k 2 + 2k(p + q) - t 2 - W 2 ) 3(k 2 - W 2 ) 3 


kg 




+ 8W 2 (k 2 - W 2 )- 


(k 2 - W 2 ) 2 (k 2 + 2k(p + q) - t 2 - W 2 ) W 2 (k 2 - W 2 )(k 2 + 2k(p + q) - t 2 - W 2 ) 


kg 


kg 


W 2 (k 2 — W 2 ) 2 4(2k§ — k 2 )(k 2 — W 2 ) 2 2k(p + q)(k 2 - W 2 ) 


kg 


kg 


+ 


kg 


The residues in the photon poles determine the value of the integration over ko variable. The integration on k can be 
performed by standard methods by means of the dimensional regularization. So in the soft approximation we obtain 
the following momentum representation of the quasipotential (52): 


V£ ft (p,q) 


(Za) 2 f 13 2 (17p 4 - 9W 4 ) 

2m 3 m 2 1 4 6 2t 2 


( 53 ) 


, 8p 6 - 4W 6 + 3p 4 (2q 2 + 3W 2 ) + 3p 2 W 2 (q 2 - 2W 2 ) , 

+ 3 ? + 

16W 8 - 29p 8 + 16p 6 (q 2 - 2W 2 ) + 3p 4 (12W 4 - q 4 ) + 4W 4 p 2 (8W 2 - 9q 2 ) , 

12t6 + 
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(44p 2 - 5W 2 ) , ^ mn , 14W 4 - 45p 4 + 38p 2 q 2 + 26W 2 p 2 , ^ mn , 
---ln(t/W) +- tto -ln(t/W)+ 


28W 6 - 8p 6 - 53W 2 p 4 + W 2 p 2 q 2 

+ It 4 


16W 4 p 2 


ln(t/W)+ 


23p 8 - 16W 8 + 16p 6 (2W 2 - q 2 ) + p 4 (9q 4 - 12W 4 ) + 4W 4 p 2 (3q 2 - 8W 2 ) 

It 8 


ln(t/W) 




The singularity of the different quasipotential terms (53) grows. The averages of the large negative powers of r < 4*- >, 
< ^ >, < > arise when we calculate the energy spectrum in the coordinate representation The transformation of 

such matrix elements was also done on the basis of relations (24)-(26), (30), (A4)-(A7) by means of the computer 
program FeynCalc for the system ’’Mathematica” |36|j37[ |. The contribution of order (Z a) 7 to the muonium (hydrogen 
atom) S-states is determined by the expression: 


AB 


soft 

n 


i 1 [605147 
+ n 9216 


(Za) V j 63565 

m 8 m27rn 3 1 288 

148257 360673 

1024 C ~~ 6144 


8675 

72 


C- 


V>( n) 


4—2 

n z 


38 o 2521 x 
y7T - y^-^(n) + 38</> (n)+ 

'3698203 7157^ 10757 ,, ' 

-1-C 4-i/>(n) 

36864 6144 12288 v v ; 


+ 



81 

64 


47 47 ,, . 

I6 C + 


where = dlnT(z)/dz, 0=0.5772156649... is the Euler constant. 


(54) 


V. DISCUSSION OF THE RESULTS. 


The calculation of the higher order perturbative corrections in a in the energy spectra of hydrogen-like systems is 
very urgent task at present due to the growth of the experimental accuracy in the measurement of the energy levels 
in many simple atomic systems. The complexity of such calculations will be enhanced with the growth of the order 
a even for one loop Feynman amplitudes. The most part of the computer programs employed for the calculation 
of Feynman diagrams is working with the mass shell particles. In the case of the bound systems the Taylor series 
of the Feynman amplitudes should be taken for the extraction of the necessary order contributions. In the present 
work we considered how to use these expansions for the calculation of the corrections (Zcc) 6 m 2 /m 2 and (Za) 7 . The 
dimensional regularization is the most important tool for the use of Taylor expansions. It gives the possibility to work 
simultaneously with the ultraviolet and infrared divergencies and to obtain Fourier transforms of the momentum 
space quasipotentials with high degree of the singularity. Many terms of the quasipotential V 2 7 obtained in this work 
in the coordinate space lead to the divergent matrix elements in the case of S-states (see Eq. (39,53)). So we have 
formulated the reqqurence relation (30) in d-dimensional space which allows to take into consideration correctly all 
coefficients in the divergent terms. Using (42) and the results of our work JlJ we can present the total contribution 
of order (Za) 6 m 2 /ni 2 for arbitrary principal quantum number n as follows (spin-dependent part of the quasipotential 
was taken into account in our calculation): 


AB[ ot 



1 

2^ 


547 

72n 3 


^(-l)”(C + *(n)-l)} 


m 2 (Za) 6 

ni2n 3 


(55) 


Numerical value of the contribution (55) to the fine structure interval 2 3 Si -4- l 3 Si of the muonium and hydrogen atom 
is equal to 0.045 MHz and 5.116 KHz. The methods for the calculation of Feynman amplitudes used in this work 
can be applied for other higher order corrections. Numerical values for the contributions of order (Za) 7 (50), (51) 
and (54) obtained in the hard and soft energy regions (5), (7) for the gross structure interval (2S -4- IS) in muonium, 
hydrogen atom and positronium are equal correspondingly 0.092 MHz, 0.130 MHz and 0.066 MHz. 
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APPENDIX A: MATRIX ELEMENTS OF THE POTENTIALS IN THE COORDINATE 

REPRESENTATION. 


To calculate the matrix elements < -B >, < > we must know the radial wave functions in d-dimensional 

coordinate space and the reqqurence relation (30). The radial wave function satisfies to the following equation: 


R + 


d— 1 


R'+ - 


1 


B 1 1(1 +d-2) 


R = 0, 


(Al) 


4 ' VA P P 2 

where dimensionless variable p = 2r\/A, B = pZa, A = —2pE. The normalized solution of this equation can be 
obtained by the standard methods. The energy spectrum is determined by (23) and corresponding radial S-state wave 
function has the form: 


Bn,i=o (r) 

where Laguerre polynomials 


2 d W d T(n) 


(2n + d — 3)r(n + d — 2) 


e 

1 . d : 


-Wr t d—2 


L d I^(2Wr), W = 


pZa 

, d—3 ' 


L n( z ) = —r z A e z -— (• 
nW nl +n V 


e z 


z A+n 


)■ 


(A2) 


(A3) 


n! dz" 

In addition to the power potential matrix elements there is need to consider the averages with the logarithmic function 
lnr: < >. Using (30) we can obtain the following matrix elements of such operators in d-dimensional space: 


lnp 

< —£■ >= - 


8e n 


(d — l)(d — 3) p 


1 4(d 2 - 4d + 5) 1 

< - > “TT 1 -“TTTT - :-7W7 < ^7 > 


(d — 3) 2 (d — l) 2 p 


4 lnp 

(d — l)(d — 3) < "p 2 "" >J 


lnp d 2 — 4d + 24 ^ 1 6 

< >= d 2 (d — 4) 2 < > + d(d-4) 


lnp 

< —t > - 


(A4) 

(AS) 


2e n (d 2 - 4d + 8) J_ 4e n 

p" 


(d — 4) 2 d 2 p 2 d(d — 4) " 2 


lnp 
< — > • 

P 2 


The averages < ^ >, < Al > are finite in d=3 space but we must find them by the use wave function (A2) taking 
into account the addenda ~ e= 4+ii in the limit d —> 3: 


lnp ?(>(n + d—2) 
< p > (2n + d — 3) ’ 


(A6) 


hip 

< —f >= 


1 


[^(d - 1) + ip(d - 2) - tp(n + d - 3)]. 


(A7) 


p^ (d — 2)(2n + d — 3) 

To note also that the matrix elements of the operators dj and (- in d-dimensional space can be calculated by means 
of the Hamiltonian operator 

d 2 (d - 1) d 


H = —- 


dr 2 


dv 




Using (23) and (A8) we obtain: 


< - > = 


2 ^ <9H ^ _ 2 a 2 

21 + d - 2 < ~dl > ~ (21 +d-2) (n +4^3)3’ 


1 


< - > u = 
r 


(n+ ^)2' 


(AS) 

(A9) 

(A10) 
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